Thermoelectric effects in strongly interacting quantum dot coupled to ferromagnetic 

leads 
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We study thermoelectric effects in Kondo correlated quantum dot coupled to ferromagnetic elec- 
trodes by calculating conductance, thermopower and thermal conductance in the Kondo regime. 
We also study the effect of the asymmetry in the coupling to the leads, which has important conse- 
quences for anti-parallel magnetization configuration. We discuss the thermoelectric figure of merit, 
tunnel magnetoresistance and violation of the Wiedemann-Franz law in this system. The results 
agree with recently measured thermopower of the quantum dot defined in a two dimensional electron 
gas. 
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I. INTRODUCTION 

Transport properties of the device consisting of the 
quantum dot attached to external leads are strongly af- 
fected by the appearance of the correlated many body 
Kondo state The phenomenon discovered long ago 
Q has manifested itself as a low temperature increase 
of the electrical resistance of diluted alloys. In quantum 
dots it shows up as an increase of conductance at low tem- 
peratures. The quantum dot devices allow the study of 
fundamental physics like Coulomb blockade phenomenon 
or Kondo effect in equilibrium and non-equilibrium con- 
ditions and in geometries not accessible in bulk systems. 
Both or one of the external leads may be normal, super- 
conducting or magnetic. In this paper we shall study the 
systems in which quantum dot is coupled with two fer- 
romagnetic leads having the same or opposite magnetic 
polarization. 

Spin polarized transport, especially single electron tun- 
neling, in magnetic nanostructures has attracted much 
interest due to its potential applications in, for example, 
spintronics and quantum computing The Kondo 
effect in quantum dots attached to normal leads (N-QD- 
N structure) has been extensively studied both experi- 
mentally and theoretically HQS- Man Y 

new ef- 
fects have been predicted and observed in the transport 
characteristics, like splitting of the zero bias resonance 
under magnetic field |5lll0|. absence of the even-odd par- 
ity effects Q or out of equilibrium Kondo effect [a, [a, El • 
The Kondo effect has also been observed in many other 
systems: single atom 0, single molecule and car- 
bon nanotubes 14]. It has also been demonstrated in 
quantum dots attached to ferromagnetic leads , where 
transport properties can, in principle, be controlled with 
aid of the electron spin degree of freedom. 

Recently we observe growing interest in electronic 



transport properties of the Kondo correlated qua ntum 
dots coupled to ferromagnetic electrodes H^|-^^. In 
such geom etry the Kondo resonance splits in parallel con- 
figuration 0, 153, 0, (however, some of the works 
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(191 122| predict no splitting), while in the anti-parallel 
configuration the Kondo effect remains virtually the same 
as for nonmagnetic electrodes. The shot noise studies 
|l8l l22j ] reveal huge differences for spin up and spin down 
electrons in the parallel alignment and no differences in 
anti-parallel configuration. 

It is well known that thermoelectric properties are the 
source of information complementary to that obtained 
from other transport characteristics ■ Thermal prop- 
erties (thermopower and thermal conductance) of the 
quantum dot coupled to the normal leads in the Kondo 
regime have recently been investigated Ther- 
mopower has been shown to be very sensitive and power- 
ful tool to study the Kondo effect. It manifests itself as 
an energy peak in the DOS slightly below Fermi energy 
and this leads to change of sign of the thermopower. 

It is the purpose of the present work to study the ther- 
moelectric properties of the quantum dot coupled to fer- 
romagnetic leads. We shall concentrate on the conduc- 
tance, thermal conductance, thermopower and related 
quantities like tunnel magnetoresistance (TMR), ther- 
moelectric figure of merit which value is direct sign of 
usefulness of the system for applications and Widemann- 
Franz ratio, which normalized value differing from 1 sig- 
nals breakdown of the Fermi liquid state. 

We show that thermopower is very different for spin 
up and spin down electrons in parallel configuration and 
is similar to the nonmagnetic case for anti-parallel one. 
However, for parallel alignment the total (spin up plus 
spin down) thermopower is very small compared to the 
anti-parallel alignment. 

The organisation of the rest of the paper is as follows. 
In Ssetion II we present the model and discuss some as- 
pects of calculational procedure. Results of calculations 
are presented and discussed in Section III. We end up 
with summary and conclusions. 
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II. MODEL AND APPROACH 

Schematic view of the quantum dot coupled to two 
leads R and L, which may be magnetically polarized 
and/or at different temperatures and voltages is shown 
in the Fig. □ 

., V 
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Vo 



FIG. 1: Schematic view of the quantum dot coupled to two 
leads R and L, which may be magnetically polarized and/or 
at different temperatures and voltages. Voltage applied to the 
junction is taken into consideration as a shift of the chemical 
potentials \il — — eV. 

We assume that the interaction energy between two 
electrons on the quantum dot is the largest energy in the 
problem and thus model the system as the U = oo single 
impurity Anderson Hamiltonian 35] in the slave boson 
representation |3t3|-[38| 
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where A = L (R) denotes left (right) lead, 



'Aktr 



(CAk<r) 



is the creation (annihilation) operator for a conduction 
electron with the wave vector k, spin a in the lead A 
and Vxkn is the hybridization matrix element between 
localized electron on the dot with the energy e a and con- 
duction electron of energy exk in the lead A. Ferromag- 
netism of the electrodes is modeled via modified conduc- 
tion electron energy exka = £Ak ± oh z where the magne- 
tization points into z direction (the same or opposite in 
both leads). 

The particle current J\ and the energy flux J ex flowing 
from the lead A to the central region can be calculated 
from the time derivative of charge and energy operator 
respectively [^J . We use relation Jq\ — J ex — HxJeX for 
the thermal flux Jq\ and express all currents in terms of 
Keldysh Green functions 01 m the standard form [2^, 



J Q x = JiJ2j ^ r ^( w )( w -^A)[G<(w) 

+2ifx(fii)lmG r a {u)] (3) 

where G r a (uj) is the Fourier transform of the retarded 
Green function (GF) G£(i,f) = i0([/ ff (t), /+(f )]+) and 
G<(w) = i(fv(t')fcAtY) is the Fourier transform of the 
lesser Keldysh GF gg. T X(T (uj) = 2irJ2 k \V X]irT \ 2 S{uj - 
fAko-) denotes the strength of the coupling between dot 
and the lead A, fx(^) = f(<-^ — fJ-x) is the Fermi distri- 
bution function in the lead A with the chemical potential 
Hx and temperature Tx- 

In general, when both the strong on-dot Coulomb in- 
teraction and the tunneling between dot and leads oc- 
cur, it is not possible to calculate Ga^ r \uj) exactly. Sev- 
eral approximation schemes have been proposed to cal- 
culate G^ t \oj). Here we use recently proposed equation 
of motion technique for nonequilibrium GF 41] . This 
technique allows calculation of both G<(ui) and G r a (uj) 
in consistent way making similar approximations in the 
decoupling procedure of both G<(lu) and G^(w). The 
approach has been successfully applied in quantum dot 
systems [33, |3^ ■ In the present case it yields 
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where T a = T La T Ra /[T La + T RtT ] 
The on-dot retarded GF reads 



and eV — /il — [i R . 
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- e a - Socr(w) - Sj(t(w) 



(6) 
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and interact- 
— . In or- 



JeX — 



du> 
2^ 



with noninteracting Eoo-(w) 

ing self-energy E Iir (w) = £ Ak —^——^^ 
der to get the splitting of the Kondo resonance in the 
presence of the ferromagnetic leads, which is consistent 
with the scaling analysis, we follow Ref. ^3] an d replace 
the bare dot energy level by e a found self-consistently 
from e a — £a + Repoo^eo-) + S/ .(e CT )] in the interact- 
ing self-energy £j CT (u;). In numerical results presented 
below we have used constant bands of width D — 100r, 
+ T/jo-) = r and use V as our energy unit in the 
following. 

Within linear response theory for the particle current 
and the heat flux one defines the conductance G as equal 
to —jrLn, thermopower is given by S = —-^fj^ and 

±k\ The 



(7) 



the thermal conductance by k = ^ — 
linear response coefficients read 



r^(u)[G<(u)+2if x (u)lmG r a {u)]{2) L n = I ^TA^lmG^uj) 
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The leads considered here are magnetically polarized. 
In the following we shall characterize the degree of po- 
larization by the parameter p = , where n a is 
the concentration of spin a electrons. The magnetization 
of both leads may point into the same direction (paral- 
lel configuration) or in opposite directions (anti-parallel 
configuration) . 
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FIG. 3: TMR as a function of the polarization factor p for a 
number of temperatures. Ed = — 1.75F, Tk(p = 0) ~ 2-10 _2 r. 



III. THE RESULTS 

In FigEl linear conductance of the system in parallel 
configuration is calculated as a function of temperature 
for a number of leads polarizations p. One observes the 
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FIG. 2: Temperature dependence of the conductance for a 
number of leads polarizations. Ed = — 1.75r. 



temperature. For small polarizations it is negative, as in 
this case Gap is not much affected by the exchange field 
while Gp is strongly suppressed due to the splitting of 
the Kondo resonance. For large values of the polariza- 
tion, TMR is positive, as in this case Gap goes to zero 
with p — > 1 while Gp tends to finite value. 

The thermal conductance k vs. temperature in the 
P configuration is displayed in the Fig. QJ As one can 
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decrease of G for larger values of p. This is due to the 
splitting of the Kondo resonance. Moreover, the decrease 
of G is not universal function of p. On the other hand, in 
anti-parallel configuration, the conductance can be ob- 
tained from that one for unpolarized leads, as it scales 
according to the relation: G{p)/G(0) = 1 - p 2 for the 
whole temperature region. 

Tunnel magnetoresistance is defined ^3 as the ratio 



TMR - 



Gp — Gap 



G 



AP 



where Gp^AP) is the conductance calculated for parallel 
(anti-parallel) configuration of the leads magnetization. 
It is shown in Fig. |31 as a function of polarization factor 
p for a number of temperatures for system characterized 
by £d — — 1.75F. Note that the Kondo temperature it- 
self is a function of p. For non-polarized leads one finds 
Tfc(p = 0) w 2 • 10~ 2 r. The dependence is not symmetric 
with respect to p = 0.5 and changes its character with 



FIG. 4: Thermal conductance of the quantum dot as a func- 
tion of temperature for a number of lead polarizations. 

see in the figure 0J thermal conductance does not change 
much with the increasing of the leads polarization P. 
This quantity is not sensitive tool for the study of the 
Kondo effect. In the case of AP configuration for any 
temperature k scales with the polarization in the same 
way as G does, namely, k(j>)/k(0) = I — p 2 . 

In the Fig. [S] linear thermopower vs. temperature is 
shown for a number of polarizations in the parallel con- 
figuration. 

The upper panel shows the total thermopower, the 
middle one the thermopower associated with spin up 
while lower one with spin down electrons. The p = 
curve corresponds to non-magnetic leads. First of all one 
can see that for p = and around T — Tk, where Tk 
is the Kondo scale, the thermopower reaches minimal 
value. It increases for elevated temperatures and eventu- 
ally changes sign. This signals the disappearance of the 
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FIG. 5: The total thermopower of the quantum dot (upper 
panel), the thermopower associated with the spin up elec- 
trons (middle panel) and that one associated with spin down 
electrons (lower panel) as a function of the temperature for a 
number of lead polarizations. 



Kondo peak. As it is well known, and in similar con- 
text has already been noted by Boese and Fazio , the 
thermopower is sensitive to the curvature of the energy 
dependence of the density of states. This curvature is 
negative at high temperatures and for Fermi level above 
the on-dot energy level, while it becomes negative at low 
temperatures, when the Kondo resonance forms slightly 
above Fermi level. 

The low temperature changes of S with polarization 
can also be understood. With increasing polarization 
in the leads thermopower decreases to zero, except at 
very high temperatures where the broad maximum (min- 
imum) is observed. Small thermopower at low T for in- 
creasing leads polarization is due to the splitting of the 
Kondo resonance by the stray fields coming from the fer- 



romagnetic leads. In this case DOS around the Fermi 
energy is almost symmetric. Moreover the thermopower 
is mostly positive in the whole range of the temperatures 
except for very small leads polarizations. It is worthwhile 
to note additional broad maximum (spin up) and min- 
imum (spin down) of the thermopower at temperatures 
below T = 0.1 for small polarizations. For such polar- 
izations the distance between Kondo resonances 6 is not 
so large (5 ~ 1) and thermal broadening significantly af- 
fects the DOS at the Fermi level, thus giving rise to the 
observed thermopower. 

For AP configuration the thermopower does not de- 
pend on the polarization and is the same as for the 
QD with non-magnetic leads. This can be easily un- 
derstood as this quantity measures the curvature of the 
DOS around the Fermi energy, which does not change 
with the polarization in this case. Only the height of the 
Kondo resonance changes. 

In the Fig. the total thermopower (upper panel), 
that in the spin up channel (middle panel), and in the 
spin down channel (lower panel) , is plotted as a function 
of the lead polarization in P configuration. Interestingly 
thermopower changes its behavior in significant way only 
for small values of p and at low temperatures it even 
changes sign. This is attributed to the splitting of the 
Kondo resonances. For small lead polarizations the split- 
ting is small but significantly influences the DOS around 
the Fermi energy. For larger p the Kondo peaks are far 
away from the Fermi energy and do not change low en- 
ergy DOS much. 

So far the results for parallel configuration of the leads 
polarizations have been presented. However, as we have 
mentioned already, for anti-parallel configuration the cal- 
culated quantities either do not depend on the value 
of the polarization and they are the same as for non- 
magnetic leads (thermopower) or can be obtained from 
the solution for non-magnetic leads due to the scaling re- 
lation 1 — p 2 (electric and thermal conductance) . This 
is easy to understand as for such polarizations the tun- 
neling of spin up and spin down electrons on the dot are 
allowed to one of the electrodes. This is true for sym- 
metric couplings. The situation is different when there is 
asymmetry in the coupling to the left and right lead. In 
this case the electric conductance behaves similarly as in 
parallel configuration (see Fig. |2J • Thermal conductance 
also does not change much with increasing of the lead 
polarization but unlike for P configuration, where it de- 
creases (increases) at low (high) temperature, it decreases 
in the whole region of the temperatures (see discussion 
below the Fig. 0}. 

On the other hand, behavior of the thermopower in 
asymmetrical AP configuration (which does not depend 
on p for symmetric couplings) is similar to the case of the 
P configuration (see Fig. |5J). In the Fig. [7| thermopower 
vs. temperature for asymmetrically coupled quantum dot 
with rLa-fTna- — 2 is plotted. At high T thermopower 
goes to zero, eventually oscillating, with decreasing tem- 
perature. 
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FIG. 7: Thermopower vs. temperature for asymmetrically 
(r_L CT /r_R CT = 2) coupled quantum dot in the anti-parallel con- 
figuration. 



tion temperature and thus one conveniently plots Z ■ T, 
which numerical value is an indicator of the systems per- 
formance, fn the Fig. (SJwe show Z ■ T as a function of 
temperature in the P configuration and note that it is 
smaller than 1, which signals limited applicability of the 
studied device. Note that polarization slightly increases 
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FIG. 6: Thermopower of the quantum dot as a function 
of the lead polarization for a number of temperatures. The 
upper panel shows the total thermopower, the middle one - 
that associated with spin up electrons, and the lower one - 
that for spin down electrons. 



Such behavior can be explained by the fact that for 
asymmetric couplings there is a splitting of the Kondo 
resonance. In symmetrically coupled quantum dot there 
is an equal number of electrons with spin up and down 
on the dot coming from different leads and one can show 
that this model can be mapped onto quantum dot with 
non-magnetic leads. Here, when TL a /Tj^ a ^ 1, due to 
different tunneling probabilities, quantum dot sees asym- 
metry in electron number with spins up and down and 
in this sense anti-parallel configuration is very similar to 
parallel one (see Fig. [SJ). 

Thermoelectric figure of merit Z — S 2 G/k is a direct 
measure of the usefulness of the material or device for 
thermoelectric power generators or cooling systems |3l[ . 
For simple systems it is inversely proportional to opera- 



FIG. 8: Temperature dependence of the thermoelectric figure 
of merit Z • T = S 2 GT/k. 

Z ■ T but it never exceeds 1. 

Again, in AP configuration the thermoelectric figure of 
merit does not depend on p, which can be easily deduced 
from the definition of this quantity. On the other hand, 
in asymmetrical AP configuration it changes little but 
unlike for P configuration, it decreases with increasing of 
the polarization. This is shown in the Fig. [5] 

Finally we discuss Wiedemann- Franz (WF) law which 
relates thermal and electrical transport via relation n = 

2 

■^TG. This law describes transport in Fermi liquid 
bulk metals and in general is not obeyed in QD systems 
where the transport occurs through small confined region 
|32l l33l| . However at very low temperatures, where the 
Kondo effect develops and the ground state of the system 
has Fermi liquid nature, the WF law is recovered. At high 
temperatures transport is dominated by sequential tun- 
neling processes leading to the larger suppression of the 
thermal transport than the electrical one. This behavior 
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FIG. 9: Temperature dependence of the thermoelectric fig- 
ure of merit Z ■ T — S 2 GT/k of the asymmetrically coupled 
quantum dot with Tlct/^Rct = 2 in anti-parallel configuration. 



is illustrated in the Fig. (|l(Jfl in P configuration (up- 
per panel) and in asymmetrical AP configuration (lower 
panel) for a number of the polarizations. 




FIG. 10: Temperature dependence of the Widemann-Franz 
ratio L/Lq = ^-(k/TG), with Lo = as a function of tem- 
perature for a number of the leads polarizations for symmetric 
coupling in P configuration (upper panel) and for asymmet- 
rically coupled dot with Flct/Yrct — 2 in anti-parallel config- 
uration (lower panel). 



In the AP configuration with symmetric couplings to 
the leads WF ratio does not depend on p due to its defi- 
nition and behavior of the quantities entering into it. 



IV. COMPARISON WITH EXPERIMENT 

It is worth to note that the thermopower of the quan- 
tum dot in the Kondo regime has recently been measured 
experimentally |43j olbeit for nonmagnetic electrodes. 
The dot containing few electrons was defined in the two 
dimensional electron gas in which also the external elec- 
trodes were defined. In the experimental situation the 
charging energies (i.e. U) were finite and varied between 
2 and 5 in units of effective coupling T. This precludes the 
direct comparison with our calculations (for p=0) as we 
have taken U = oo limit. Experimentally one applies the 
bias Vsd between source and drain and measures ther- 
mopower as function of gate voltage V g . Theoretically 
this dependence can be modeled by plotting S as func- 
tion of on-dot energy level Ed- Experimental data, shown 
in figure 4 of reference 43] qualitatively agree with our 
calculations for p = presented in figure ITT1 



IT, 0.0 




FIG. 11: The linear response thermopower of the quantum 
dot coupled to non-magnetic leads as a function of the po- 
sition of the on-dot electron energy calculated for different 
temperatures (in units of T): T = 10" 3 , 10" 2 , 5 • 1CT 2 , 10" 1 
and 2 ■ HT 1 . 

To understand the results, as the first approximation 
one can use so called Mott formula, which states that 
thermopower is proportional to the logarithmic deriva- 
tive of the conductance with respect to the energy evalu- 
ated at the actual Fermi energy. The Kondo effect shows 
up as additional maximum appearing in conductance via 
quantum dot when the temperature decreases. It is es- 
sentially due to Abrikosov-Suhl resonance in the density 
of states, which appears at low temperatures and is lo- 
cated slightly above the Fermi level of the leads. The 
slope of conductance thus changes and the thermopower 
changes sign. The similarit y b etween theoretical Fiallll 
and experimental Fig. 4 of |43j results is very encourag- 
ing. However for qualitative comparison one has to take 
few energy levels on the dot, assume finite U values and 
selconsistently calculate the shift of Ed with changing gate 
voltage and nonlinear conductance and thermopower for 
actual value of V g and Vsd- This is outside the scope 
of the present paper. The detailed comparison between 
experimental data and calculations will be the subject of 
future work. 
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V. SUMMARY AND CONCLUSIONS 

In summary, we have studied thermal properties of the 
strongly correlated quantum dot coupled to the ferromag- 
netic leads in the Kondo regime. We have found that 
thermopower is strongly suppressed at low temperatures 
due to the splitting of the Kondo resonance in parallel 
configuration of the lead polarization. In anti-parallel 
configuration Kondo effect behaves in the way similar to 
the system with non-magnetic electrodes so the results 
do not depend on the value of the polarization in the 
leads. Moreover we have shown that asymmetry in the 
coupling to the leads in anti-parallel configuration has 
important consequences as it lifts spin degeneracy on the 



dot thus leading to the suppression of the thermopower at 
low temperature, similarly as in the parallel polarization 
configuration. Finally we have checked the Wiedemann- 
Franz relation which does not hold in general, but simi- 
larly as for QD with non-magnetic leads it is recovered at 
low temperatures where the Kondo effect develops. The 
results qualitatively agree with experimental data. 
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